On the signature of tensile blobs in the scattering function of a stretched polymer 
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We present Monte Carlo data for a linear chain with excluded volume subjected to a uniform 
stretching. Simulation of long chains (up to 6000 beads) at high stretching allows us to observe the 
signature of tensile blobs as a crossover in the scaling behavior of the chain scattering function for 
wave vectors perpendicular to stretching. These results and corresponding ones in the stretching 
direction allow us to verify for the first time Pincus prediction on scaling inside blobs. Outside 
blobs, the scattering function is well described by the Debye function for a stretched ideal chain. 
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The notion of blobs in polymer physics is a widespread concept which indicates a change of fractal dimensionality 
of the polymer coil across a characteristic length £ (the blob " size" ) intermediate between the radius of gyration of the 
chain and the local Kuhn segment length. Thermic blobs in highly dilute solutions and swollen blobs in semi- dilute 
solutions have been observed by careful SANS experiments (!]] focused on the asymptotic regime of the single chain 
structure factor S(q) oc q~ x l v where v is the scaling exponent. On semi-dilute solutions in good solvent, Farnoux 
et al. |l) observed a clear change from the excluded volume statistics (y — 0.6) at length scales characteristic of the 
interior of the blob to the ideal chain statistics (y = 0.5) over larger length scales where intramolecular excluded 
volume forces are screened by the presence of other chains. The size f w 1/q* of the "swollen" blob is naturally fixed 
by the cross-over wave vector q* separating the two regimes. 

Dilute solutions of polymers under flow and stretched chains are other notorious examples in which the blob concept 
has been invoked under the terminology of "shear" blobs and "tensile" blobs respectively. However SANS experiments 
are more difficult to perform in those cases partly because effects are anisotropic and S(q) thus depends upon the 
direction of q with respect to the field. 

Onuki has conjectured that shear blobs could develop in long linear chains of N bonds in dilute solutions subjected 
to high shear rates ||] . The shear blob size would correspond to the size of a polymer segment of length n (n < < N) 
such that the longest relaxation time of this segment is of the order of the inverse shear rate. The mere existence of 
such blobs has not yet been tested experimentally although the scattering function of polymer solutions under shear 
has already been probed by SANS experiments in the power law regime H . Detecting blobs would require measuring 
the single chain structure factor of very long and highly stretched chains, a measurement made difficult by a very 
unfavorable signal to noise ratio at the very low concentrations required to get by extrapolating "infinite dilution" 
results (§. 

In this letter, we report Monte Carlo calculations on the structure factor of very long chains (up to 6000 monomers) 
stretched at their ends by equal but opposite forces f which provide direct evidence in favor of the tensile blob concept 
introduced 20 years ago by Pincus and de Gennes H|| . These authors suggested that when polymers in good solvent 
are stretched by sufficiently high tensile forces, the coil behaves like a stretched ideal chain formed of tensile blobs as 
basic units. Inside the blobs, good solvent statistics persists. This picture rests on the assumption that the strong 
elongation of the chain in the field direction prevents monomers distant along the polymer contour to interact with 
each other directly. In this picture the blob size £p depends on the force intensity as £r = kgT / f where fee is the 
Boltzmann constant and T the temperature. 

Macroscopically the onset of tensile blobs above a threshold value of the reduced force 77 = R /£t (Ro is the 
unperturbed chain end-to-end distance) results in a change of the elastic behaviour of the chain. For small stretching 
the elongation Rf of the chain, defined as the average projection of the end-to-end vector in the field direction, follows 
a linear Hooke's law Rf = R rj/3. In the same 77 regime, fluctuations of the square end-to-end distance parallel and 
perpendicular to the force remain at the equilibrium value. Above a certain value rj* one obtains the following power 
laws 

R f = RoBrf' 3 (I) 
5R\ = i^C|,rr 1/3 (2) 
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Rl = R 2 C ±V -^ (3) 

where B, C\\ and C± are model dependent constants (independent of N and 77) and the exponent values 2/3 and —1/3 

'I 



arise from the expected linear dependence of Rf, R\ and SR?, on N at fixed / in this regime. In what follows, we 



define R\ as the fluctuations in any direction perpendicular to stretching. 

A direct test of the elastic behavior, as given by eqs. (|l|, ||) came long ago from the precursory MC simulation of 
Webman et al. which observed the threshold to be located around r/* — 2. More recent MC study || has been 
focused on the cross- over region in the attempt to test renormalization group theory predictions || . 

We concentrate on the high stretching regime and we provide for the first time numerical data for the single chain 
structure factor. Chain lengths of many thousands beads allow us to provide the first experimental test of Pincus 
scaling J5| . Data for q perpendicular to the force exhibit a well defined crossover between ideal and self-avoiding chain 
statistics in a way analogous to the semidilute case studied by Farnoux et al. jij. In the direction parallel to the 
field and for large enough q's, the scattering function follows the usual unperturbed chain asymptotic behavior when 
plotted versus an effective scattering vector taking into account the change in metric due to stretching, as predicted 
by Pincus Q. In the low q regime, the generalized Debye function for a continuous gaussian chain [|l0| very accurately 
explains all the spectacular features of the scattering function (see figure 4) . 

Our Monte-Carlo study is based on a polymer model which is a simple necklace of hard spheres connected by 
rigid bonds taken as length unit. The model is characterized by a hard sphere diameter of a = 0.65, a value which 
ensures a fully developed excluded volume statistics |ll| . Equilibrium calculations for this model lead to the scaling 
law: R = 1.011 TV ' 60 . We combined a reptation algorithm with the Configurational Bias Monte Carlo (CBMC) 
method jl2| suitably adapted to include the external field influence. A new conformation is attempted by suppressing 
a segment of k beads from one end of the chain chosen randomly and constructing a new segment of same lenght on 
the opposite side. Leaving out details about the evaluation of the Roscnbluth weights appearing in the new segment 
acceptance probability |l2|], let us just mention here that the set of p trial directions for the single bond reconstruction 
has been extracted from a probability oc exp{f coscj) / ksT) where <p 1S the angle between the trial direction and the 
external force. The number of such trial directions has been chosen in a way to minimize the CPU time between 
two independent configurations. Two configurations generated during the MC Markov Chain are supposed to be 
independent when they are separated by (N/2k) 2 successful reptation moves of k beads. Our results are based on 
statistics of 4000 independent configurations for the chains of length N < 6000 and w 1200 independent configurations 
for the N = 6000 case. The optimal number of trial directions p is found to increase slightly with the chain lenght 
(typical value is p = 5). Chain lengths of N — 512, 1024,2048 and 6000 bonds were studied at various stretching 
forces in the range 2 < 77 < 40.36. 

Figure 1 shows a Kratky plot of the structure factor of a 6000 bead stretched chain, namely q 2 S(q) versus q for the 
scattering vector pointing perpendicularly to the stretching direction. For the three particular values of 77 considered, 
the crossover from ideal to excluded volume behavior clearly takes place within the q range corresponding to the 
universal regime of S(q) for unstretched chains. The latter is known to cover a range 3/i? g < q < A -1 where R g is the 
equilibrium radius of gyration and A the Kuhn segment length ||l3|| . In the present case such range cover two decades 
from q m i n ~ 0.04 up to q max ~ 4.0. In figure 1, we also report the Debye function in the perpendicular direction |ic| l 

S{qx) = ^££{e-*--l + X x ) (4) 

where Xj_ = q 2 L R 2 L /2 and values of R\ are directly obtained in the simulation. The Debye function is seen to follow 
very closely the chain structure factor up to the crossover region beyond which all three experimental curves converge 
towards a unique power law behavior in q 1 / 3 (represented in the figure as a straight line) which is reminiscent of 
the equilibrium universality S(q) oc q~ 5 / 3 . The observed asymptotic cross-over wave vectors q + (rj) (taken as the 
intersection value between the excluded volume power law and the Debye functions) should be related to the blob 
size £t- Indeed one can easily verify that q + {ri) — 4.5/£r(?7) ^ 27t/£t(?7)- 

An important issue is about the universality of S^q, f , N) . In figure 2, we have plotted the form factors S(q) / (N + 1) 
for various chain lengths at the same reduced force r; = 7 versus qRo- In each specific direction we get a dis- 
tinct universal function, which suggests an extension of the well known equilibrium universality of S(q) towards 
S(q±) — (N + l)f±(q±Ro,r]) and S(q\\) — (N + l)h(quRo,r}). Instead if we compare S(q) data for different chain 
lengths at the same absolute force intensity /, all curves are found to merge for q > q* as it should be. 

Along the lines of previous calculations, we have analyzed the end-to-end vector statistics of stretched chains 
and we confirm that, in the high stretching regime, the elongation as well as the squared fluctuations follow the 
universal laws eqs.(Q, g, ||) in terms of 77. The model dependent parameters are B = 0.46, C\\ = 0.31 and C± = 0.46. 
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Figure 3 shows the systematic evolution of S(q±) and S(q\\) as rj increases for N — 6000. In the direction of the 
force we observe the appearance of a shoulder which becomes more and more pronounced as the chain extends. As r) 
gets even larger, oscillations develop beyond the initial decay. The origin of this behavior can easily be understood 
by comparison with the structure factor for a rigid rod of length L with N + 1 scatterers uniformly distributed which 
is, for q parallel to the rod 0, 

S(q) = 2{ *+ 2 1 \ l-cos(qL)) (5) 



In figure 3, the comparison with the rigid rod is made for r) — 40.36, using L — R Br] 2 / 3 (see eq.(jlj)). At small q 
the chain can be described as a rod of fluctuating length, a point not totally surprising if we note that for N = 6000 
and rj — 40.36, the average extension is one sixth of the polymer contour length! The rigid rod model becomes 
obviously more and more inaccurate for increasing q, that is on the length scales where fluctuations start dominating 
the systematic effects. Further, we note that beyond q = 10/£r, indicated by vertical bars in the figure, the structure 
factor remains isotropic. 

Figure 4 shows the structure factor for gn in the N — 6000 case for four values of rj = 0.00, 7.00, 13.34, 40.36. Data 
for rj ^ are also reported versus the effective wave vector q defined as 

q 2 = q 2 + w cos 2 6 (6) 
vr 

where 9 is the angle between the observation direction and the external field p). For cosO ^ 0, q is a non linear 
combination of the bare q vector and a measure of the systematic stretching effect. For q 3> ^y 1 , q — q and we recover 
the typical power law of the equilibrium scattering function. In an intermediate region q** < q < 10/£t the curve 
bending of S(q\\) produced by the external force and directly related to the presence of a finite systematic stretching, 
is compensated by plotting S(q» ) versus q. This brings back the data on the equilibrium power law p] 

S(q\\)~q- 1/U (7) 

Finally in the region q < q** , q becomes almost constant so that the data exhibit an apparent divergence when plotted 
against q. 

As for the perpendicular direction, the continuous gaussian chain model can be exploited to explain the S(q\\) at 
low q. The generalized Debye function is [|T(| 

/ er x « - 1 + X\\\ 

S(q\\ ) = 2(N+ 1)SR - 2 1 (8) 



where Xn is a complex variable defined as 



X^^Cfr^ + ^RoBr?'* (9) 

In figure 4 we report the Debye behavior (dashed lines) for the three values of rj and for values of the elongation and 
the parallel fluctuation given by the observed power laws behaviors of eqs.(QJ^). The agreement between the ideal 
chain prediction and the experimental data is amazingly good up to q values close to q** . 

In conclusion we have simulated long linear chains stretched by a uniform force. We have studied the strong 
stretching regime where tensile blobs dominate the elastic behavior of the chains. We have performed a careful 
analysis of the chain scattering function and thanks to the large dimension of our chains we have been able to 
directly observe the blob signature as a cross- over in the universal behavior of the scattering function in the direction 
orthogonal to the force. We have provided the first "experimental" evidence that Pincus scaling in terms of an effective 
wave vector || describes accurately the entire excluded volume region of the scattering function (inside the blobs). 
At larger length scales the generalized Debye function beautifully explains the scattering data, in particular at very 
high extension when the chain behaves as a collection of blobs aligned as a linear rod. The picture of an ideal chain of 
Nb blobs of size £t stretched by a force / is however only qualitatively supported by our data. Such a model requires 
R± = SR 2 = (l/3)Nb£r together with Rf = (l/3)A"b^(//fcsT) which according to the observed behaviour of eqs.([l], 
||, |3|) demands C± — C\\ = B. In our results, C± — B is indeed observed but C\\ ~ 2/3C±. 

The relevance of our results with respect to the experimental situation is threefold: first, it should be noted that 
single chain stress-strain curves are nowadays directly measured in beautiful experiments on single DNA molecules p4[ . 
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However, similar experiments on flexible uncharged single polymers in good solvent still appear unfeasible. The second 
and more traditional experimental source of information on stretched chains would be scattering measurements on 
highly diluted solvated gels. However, well known network heterogeneities of such samples will transmit stress on 
junction points in a uncontroled way, making it very difficult to relate observed spectra to single chain statistics. 
Finally, it appears as most likely that the first experimental test of the features of the chain structure factor discussed 
in the present letter will be provided by SANS experiments on sheared polymer solutions. 

We have the pleasure to thank G. Destree for unvaluablc technical assistance. We thank P. Lindner, D. Richtcr 
and J. Titantah for illuminating discussions. We benefited from a generous INFM computer grant on the CRAY-T3D 
parallel computer at Cineca, Bologna (Italy). 
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FIG. 1. N=6000. Kratky plot of S(q) q 2 versus q for r\ — 7 (circles), 13.34 (squares), 40.36 (crosses). Dashed curves are the 
predictions of the Debye function and the straight line represent the excluded volume behavior S(q) q 2 ~ q 1 ^ 3 . 

FIG. 2. Universal behavior of S(q)/(N + 1) versus qR at r\ = 7. Various chain lengths are shown for both q parallel (lower 
curve) and perpendicular (upper curve) to the external field. N=512 (circles), 1024 (squares), 2048 (crosses) and 6000 (stars). 
Equilibrium behavior for N=1024 is given by the thick line. 

FIG. 3. N=6000. S(q) for q perpendicular (upper curves) and parallel (lower curves) to the field for rj — 7 (circles), 13.34 
(squares), 40.36 (crosses) versus the bare wave vector q. The continuous line is the results of the rigid rod model mentioned in 
the text. 

FIG. 4. N=6000. S(qu) for rj = 7 (circles), 13.34 (squares), 40.36 (crosses) versus the bare wave vector q (lower curves) and 
the rescaled wave vector q (upper curves). For each value of ij, the lower boundary q**(rj) for the validity of Pincus' scaling is 
indicated by a vertical bar. Dashed curves are the predictions of the Debye function. Note that they are very accurate up to 
q**(r]) which is therefore the crossover from ideal to excluded volume behavior. 



4 



